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SUMMARY 


The frequency spectrum at 0°K and the specific heat at low temperatures have been calcu- 
lated in three different approximations. The Born model is shown to give results which are defi- 
nitely not in accordance with experiments. A second approach which has some theoretical 
support in the Heitler-London treatment of ionic crystals gives considerably better results. 
Finally, this approach has been supplemented by taking the polarisability of the ions into account. 
The polarisation is shown to have a strong influence upon the frequency distribution. The results 
in this approximation are those which best agree with the observed specific heats, although some 
discrepancy still remains. 


Introduction 


The frequency distribution of thermal vibrations and the specific heat at low 
temperatures for KCl and NaCl have earlier been investigated by Kellermann [7] 
and Iona [6]. Both investigations were based on the Born model, i.e. a long range 
Coulomb attraction between the ions and a short range repulsive central force. The 
: agreement between the calculated and observed specific heat was remarkably good 
in both cases. This fact has often been considered as a confirmation of the approxi- 
mate validity of the assumptions inherent in the Born model. The significance of 
their results must, however, be questioned for the following reason. The force constants 
associated with the repulsive forces were in both investigations chosen so as to fit 
the elastic constants at room temperature. The calculated frequency distribution 
must therefore be considered as an approximation to that at room temperature. The 
comparison with experiments, however, concerns the specific heat at low tempera- 
tures (0° — ~100°K). Now, if the force constants are instead chosen so as to fit 
‘the elastic data at 0°K, the frequency of each individual mode of vibration will 
change significantly, and it may be expected that the corresponding frequency 
distribution, which would be an approximation to that at low temperatures, will 
bear no close relation to that calculated in the earlier work. This has here been con- 
firmed by actual calculation and the results given later in this paper show that the 
results obtained are not in agreement with experiments. 
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A natural starting point for an investigation of the lattice vibrations would be 
provided, if the many-electron problem for an ionic solid had been solved with 
sufficient accuracy. As this is not the case one is forced to proceed in a phenomeno- 
logical manner and one is then naturally led to consider the information brought 
out in several investigations of the long polar oscillations of ionic solids,! which can 
be discussed using macroscopic theory. In first place it has been found that the elec- 
tronic polarization is of essential importance; the main effect being to cut down con- 
siderably the limiting frequency of the longitudinal vibrations. It has been anticipated 
in the literature that the polarisation may effect also other modes of vibrations in 
a singificant way, although no actual calculations of the effect have been done up 
to now. Unfortunately, we have not been able in this investigation to include a 
satisfactory treatment of the polarisation, but have had to limit ourselves to the 
rather idealized case, where the lattice is considered as an assembly of point-dipoles. 
We believe that such a treatment cannot be very accurate, but hope that it will at 
least give the correct trend. 

A second effect of importance is concerned with the effective charge. This concept 
was first introduced by Lyddane Sachs and Teller [10] and was later investigated in 
more detail by Szigeti [14], who derived by macroscopic considerations a formula 
expressing the effective charge in wellknown empirical data. He also pointed out 
that deviations from the ionic charge might be due to a “‘distortion dipole” caused 
by the mutual distortion of neighbouring ions. The study of the limiting vibrations 
gives, however, no idea how the effect should be treated for a general mode of vibra- 
tion and to proceed it therefore seems necessary to consider a specific model for 
the interaction between the ions. In a study of the Heitler-London approximation 
in connexion with lattice vibrations by one of us? it was found that the concept of 
effective charge enters without any ad hoc assumptions and agrees qualitatively 
with the results of Szigeti. Quantitatively, however, the accuracy of the Heitler- 
London method is not sufficient to make the method useful for calculations, and 
therefore the treatment was simplified so that the most important parameters 
could be determined from empirical data instead of calculated using electronic wave 
functions. For this investigation, however, we have preferred not to make full use 
of the formulas in ref. [9], but instead to use a simplified version of that scheme, 
which will be briefly described in the next section. 


Brief description of the theoretical basis of the calculation 


The frequencies of each mode of vibration are found from the plane-wave solutions 
of the equation 


ay = ie al: (1) 


Here the index | refers to a unit cell, & distinguishes between the particles within 
a cell, m,, is the mass of the particle at the position &, and uj, denotes the displacement 
of the particle (1, &) from its equilibrium position. ®,,: is a tensor depending on the 
interaction between the ions and is found from the lattice potential. As the contribu- 
tions to ®}, due to effective charge modifications and polarisation are discussed in 


1 We refer to [3] for a comprehensive review of this subject. 
2 §. O. Lundqvist [8], [9]. 
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detail in ref. [9], we shall only briefly mention the further simplifications introduced 
here. The simplification is entirely concerned with the contributions associated with 
the effective charge. According to the Heitler-London scheme the effective charge 
is described as a many-body effect, giving rise to terms in ®;,): of a rather complicated 
appearance. In order to simplify the calculations we have neglected most of these 
terms and only retained the simplest one, having the form of a Coulomb interaction. 
In this approximation ®;;- can be expressed as a sum of three parts 


Dj = OD + PD ie + PBL, (2) 
which will be discussed separately below. 
The first part with components 
lv’ a 1 
COD KK’ = 7 &, Ey € : ; ( ; ) (3) 
2 pa oat | 
m B 


represents the modified Coulomb interaction. Here ¢, and ¢,, denote the valencies 
of the ions k and k’. ¢ is the elementary charge and the indices « and f refer to the 
Cartesian components. The parameter 7 is determined from the average of the elastic 
constants ¢,, and ¢,, according to the formula! 


_ (Get C4q) a4 
dA 0.69540 a 


where a is the interionic distance. It should be noted that the deviation from the 
classic Coulomb interaction in (3) is due to a many-body effect in the interaction 
potential and does not correspond to a modified two-body Coulomb interaction. We 
also point out that 7 is an approximation to (e*/e)?, e* being the effective charge 
appearing in the polarisation term (8) discussed later. In general this approximation 
is not satisfactory but for KCl 7 agrees fairly well with the estimate given by Szigeti, 
and for NaCl the use of (3) anyhow gives an improvement. 

The second part 

Ww a ; 
eae Sap yp {V (|r — Pie |)} (5) 

O ackk’ 0 ackk’ 

% B 
is mainly due to the repulsive interaction between nearest neighbours.? There would 
of course be no difficulty to include also a central interaction between next nearest 
neighbours except for the fact that the number of empirical constants is insufficient 
to determine more parameters. For the calculations we write (5) in the form 


il’ i. LP 
PPD KK) = wkk! xkk’ Y + bap P. (6) 
ap a B 


@ and P are determined from the equilibrium condition and the elastic constant 
¢,, according to 
1 Because of the approximation introduced here the Cauchy relation c,, = ¢4, is automatically 


fulfilled, so in all calculations we replace either of them by their arithmetic mean. 
2 According to ref. [7] and [8] V(|ri, =r) also contains a part derived from a three-body 


potential, but which here acts as an ordinary central potential. 
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e2 
ey =aQ—7 7. 1.569 (7) 
and ‘ 
ne 
Pin 
a3 -0.2913 7) 


Before giving the expression for the third part in (2), which describes the effect 
of the polarization, we introduce two matrices T and @ having the elements 


iV’ Ge | on 
[kk = ; 7 : 
Ul il |r! _y! a 
O akk’ Oackk’ MOR OR 
CD a8 


and 


ll’ 
akk! = ot Or Onn: O 28 
ap 


where «, denotes the polarisability of the ion k. 
In this notation 


pol ae = 42 $7 = —j LT 
Dkk’ = ee, e*" {T (a > —T) * Thee (8) 
ap op 


where e* denotes the effective charge. These coefficients are not easy to compute 
according to (8), but what is needed for the calculations is actually the associated 
Fourier matrices, and for these the problem reduces to operations on 6 x 6 matrices 
only, which can be handled by standard methods. 

The calculations have been performed according to the following three approxi- 
mations: 


(a) The Born model in which the polarisation term (8) is neglected and the para- 
meter 7 in (3) equals one. 

(b) An intermediate approximation in which the polarisation is neglected but the 
parameter 7 is determined from (4). 

(c) Both the polarisation and the modified Coulomb interaction are taken into 
account according to (3), (4) and (8). 


The calculations were done for 89 non-equivalent wave vectors evenly distributed 
in 1/48 of the reciprocal lattice cell, giving totally 534 frequencies. Tables were 
prepared for the contributions from the Coulomb and repulsive forces, and the cal- 
culation of the polarisation term and the final solution of the secular equations were 
done using the electronic computer BESK in Stockholm. Approximate frequency 
distributions {(@) have been determined by counting the number of frequencies in 
intervals of the length Aw = 0.2-10- radians/sec. The distribution functions were 
normalized such that {f(w)d@ = 6. Because of the rather small numbers of frequencies, 
we have preferred to consider the histograms, instead of drawing continuous curves. 
Due to the same limitations we have not found it possible to make a closer study 
of the nature of the singularities in the distribution functions. 

The specific heat was calculated from the formula 


we 
Cr=kN| ppl (odo 7 
“Ya (9) 
T° kT 
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where N is the number of the oscillators, # Bolzmann’s constant, on =h/2a where 
h is Planck’s constant, w are the frequencies of the oscillators, and /(q) is the distribu- 
tion function of lattice vibrations. Part of the uncertainty due to the approximative 
frequency distribution functions could possibly have been eliminated by expressing 
Cy as an integral over the reciprocal lattice cell, but due to the tentative nature of 
our calculation, we preferred to use (9), which leads quicker to results. 


197 


S. O. LUNDQVIST et al., Lattice vibrations of KCl and NaCl at low temperatures 


The frequency distribution and specific heat of KCI at low temperatures 


In this calculation the following set of empirical constants have been used 


KCl 
dy (interionic distance): 3.115- 10-8 em 
Cra: 4.95 - 10" dyn/cm? 
Cio: 0.6: 10" dyn/cm? 
Cag: 0.669 - 1044 dyn/cm? 
Ok: 1.13 - 10-74 em3 
aon! 2.92 - 10-24 cm 


The elastic constants have been taken from Durand [4] and the polarizabilities 
from Shockley et al. [11]. The parameter 7 was found to be 0.74. We note that 73 
is an approximation to e*/e. Szigeti finds e*/e =0.80 which is rather close to the 
value 7? = 0.86. As Szigetis values is found from empirical constants at room tempe- 
rature the figures are not quite comparable, and we preferred to use the value e*/¢ = 
0.86 in our calculations. 

The frequency distributions are given in Figs. 1a, 1b and 1c. A comparison between 
the curves shows that one important effect of the effective charge and polarisation 
is to shift the high-frequency part towards lower frequencies. If Ib and Ic is compared 
one finds that not only the high frequency part of the spectrum is affected by the 
polarisation but also the intermediate frequency part (~2-10-!% rad./sec.). In order 
to see more directly how the frequencies of the individual modes are changed in the 
different approximations we have given the curves wm =w(k) for waves along a 
crystal axes in Fig. 2a, b. The frequencies according to approximations b and ¢ are 
given in Table 1. 

The results for the specific heat are given in Fig. 3, where the Debye temperature 
(py is plotted against the absolute temperature 7’. The results obtained from the 
Born model are, as seen from Fig. 3, far off from the experimental results.1 Already 
the inclusion of a modified Coulomb interaction changes the result significantly, 
but gives still too high values of 6p over the whole temperature interval. The effect 
of the polarisations is roughly to shift the whole curve downwards, but it seems 
clear that our treatment of the polarization overestimates the effect. 


The frequency distribution and specific heat of NaCl at low temperatures 


The empirical constants used in our calculations for NaCl are 


Gg weeloy LOz? em 

C4; 0.75 +10" dyn/em? 
Cy:  0.986-1011 dyn/cm? 
Cy: 1.327-10" dyn/em? 
Ona: 0.28 -10-%4 cm? 
Og: » 2:92) -10-%4 em? 


1 Note that our treatment differs somewhat from Iona’s, which also includes a central interac- 
tion between next nearest neighbours. 
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Fig. 2a. The transversal branches of vibration Fig. 2b. The longitudinal branches of vibration 
for KCl a) -—- b) —— ec) for KCl a) b) c) 


) ie ae 


The Debye characteristic temperature for KCl. 
Fig. 3. a) — b) ec) . The experimental points are taken from Berg and Morrison [2]. 
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The values of the elastic constants are those given by Overton and Swim [10] and 

the polarisabilities are taken from Shockley et alien [11]. The parameter 7 was 

found to be 0.87. This gives a rather poor estimate of the effective charge, 7! = 0.93, 

so for e*/e we have used the value 0.74 found by Szigeti [12]. 

_ The frequency distributions are reproduced in Fig. 4a, b, c. The changes in 
the distribution function in the different approximations are qualitatively the same 
but are more accentuated for NaCl. Especially notable is the strong suppression of 

the high frequency part in (c) compared with (a) as well as the shift of the position 

and height of the peak occurring approximately at the residual ray frequency. The 

relation w =w(k) is given for waves along a crystal axes in Fig. 5, a, b, and shows a 
significant change due to the polarisation for all modes of vibrations except the 

acuostical vibrations of long wavelength. The frequencies according to approxima- 

tions b and ¢ are given in Table 2. 
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Fig. 4a. The frequency distribution function for 
NaCl according to approximation a. 
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Fig. 4b. The frequency distribution function for 
NaCl according to approximation b. 


Fig. 4c. The frequency distribution function for 
NaCl according to approximation c. 


The curves of 0p versus T is shown in Fig. 6. As for KCl the Born approximation 
is in bad agreement with experiments; the inclusion of a modified Coulomb interac- 
tion alone gives some improvement, and together with the polarisation one obtains a 
strong shift of the whole curve, but the result indicates that the effect of the polarisa- 
tion is somewhat overestimated in our treatment. 
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Fig. 5a. The transversal branches of vibration Fig.5b. The longitudinal branches of vibration 
for NaCl a) —b) c) for NaCl a) b) c) 
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The Debye characteristic temperature for NaCl. 
Fig. 6. a) -—— b) c) . The experimental points are taken from Clusius, Goldmann, 
Perlick [4] and Martin [11]. 


Conclusion 


The results show that for both crystals the qualitative features are the same in all 
three approximations. For both KCl and NaCl the Born model gives results which 
are not in agreement with experiments. The approximation where only the Coulomb 
interaction is modified as well as that where the polarisation is included agrees better 
with experiments; although in none of the cases the agreement is good. This is cer- 
tainly due to the very approximate theory on which the calculations are based. 

“Nevertheless the results show clearly that both the effective charge modification and 
the polarisation are of importance and have to be considered in calculations of the 
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vibrations of ionic solids. It also seems rather satisfactory that the experimental 
eurves for both crystals are bracketed in by two fairly simple approximations, and 
that, in spite of the approximate nature of the treatment, the accuracy of the results 
in both approximations seem to compare well with that obtained in several similar 
investigations of non-ionic solids. Work is now under progress in order to find a 
more satisfactory theoretical basis for further investigations. 


We wish to thank the mathematical staff at the electronic computer BESK in Stockholm for 
their helpful cooperation. The financial support of he Swedish Natural Research Council is grate- 
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‘ Appendix 
able J. Values of the frequencies for KCl. 


For each point the upper value belongs to the approximation b and the lower value to the ap- 
proximation c. pz are the components of the wave vector in units of }- (1/2a,). 


2 2 2 2 | s = s 8 . = 2 
Ee mb e|rer: | & ee le oe 
- a = ‘2 = ial a ° ial s 

3 3 3 8 8 3 3 8 8 8 S) 8 


00} 4.50] 2.95] 2.95} 0.00} 0.00] 0.00 SDs Ze! 3:69.) 2-16) *2.04,|" 2.83 | *1,64] 1.33 
4.18} 2.68] 2.68] 0.00] 0.00] 0.00 3.45 | 2.55] 2.38| 2.25] 1.56] 1.33 
00 | 4.44] 2.95] 2.95] 0.61] 0.22] 0.22 430] 3.57| 3.02] 2.39] 2.30] 1.80} 1.24 
4.11} 2.68) 2.68] 0.61] 0.22] 0.22 3.31} 2.69| 2.30] 2.25] 1.75] 1.06 
10] 4.43} 2.95] 2.90] 0.72] 0.58] 0.32 43 1| 3.59] 2.91] 2.39) 2.33]. 1.83] 1.32 
4.11) 2.68 | 2.66) 0.72] 0.57) 0.31 3.35] 2.63] 2.30] 2.29] 1.77.| 1.20 
11| 4.41) 2.90] 2.90] 0.82] 0.61] 0.61 43 2] 3.58] 2.70] 2.46] 2.36] 1.89] 1.52 
4.11} 2.65] 2.65} 0.82] 0.61] 0.61 3.38| 2.50} 2.39] 2.28] 1.82] 1.48 
00] 4.28] 2.96] 2.96] 1.21] 0.44] 0.44 433] 3.48| 2.64] 2.50} 2.32] 1.98] 1.74 
3.92] 2.68] 2.68) 1.19} 0.42] 0.42 3.32| 2.58] 2.36] 2.23] 1.91] 1.74 
10; 4.28] 2.96] 2.87] 1.24] 0.73] 0.50 440] 3.34] 3.03] 2.43) 2.21] 2.09] 1.45 
3.96| 2.68} 2.64] 1.22} 0.72] 0.49 3.11} ° 2.71 |* 2-42 | °2.17-] 1.99] 1.19 
yl 4:28) 2.91) 2.85) 1:27) 0.81] 0:71 441) 3.38] 2.91] 2.49] 2.21] 2.09] 1.52 
3.98] 2.65) 2.62] 1.26] 0.80] 0.71 3.16] 2.62] 2.46] 2.16} 2.00] 1.31 
20) 4.20| 2.97) 2.74] 1.39] 1.13] 0.66 442] 3.38] 2.68] 2.63} 2.20] 2.09] 1.69 
3.90} 2.68] 2.57] 1.39} 1.12] 0.63 3.20} 2.59| 2.45] 2.15] 2.02] 1.59 
21) 4.20| 2.89} 2.74] 1.45] 1.15] 0.85 443] 3.29] 2.85] 2.41] 2.20] 2.09] 1.90 
3.92| 2.64] 2.57] 1.44] 1.13] 0.84 3.15] 2.78] 2.28] 2.13] 2.05] 1.87 
22 | 4.13) 2.75| 2.75] 1.63) 1.19] 1.19 444] 3.17] 3.02] 2.20] 2.20] 2.09] 2.09 
3.88 | 2.56] 2.56} 1.62} 1.19] 1.19 8.03 | ~ 2°96) 2.11 |) 2.11 | 2.07) 2.07 
0 0} 4.01} 2.97] 2.97] 1.75] 0.64] 0.64 5 00} 3.31] 3.00] 3.00} 2.60] 0.97] 0.97 
3.62] 2.68] 2.68] 1.70] 0.60] 0.60 2.84) 2.68] 2.68} 2.31] 0.85] 0.85 
10] 4.03] 2.98] 2.86] 1.76] 0.86} 0.70 5 10] 3.36] 3.01] 2.86] 2.60] 1.13] 1.03 
3.67] 2.68] 2.62] 1.72] 0.85] 0.65 2.97} 2.68| 2.55] 2.31] 1.05] 0.90 
11] 4.04] 2.92] 2.83] 1.78] 0.96] 0.84 5 11] 3.40} 2.94] 2.80] 2.60] 1.26] 1.10 
3.71] 2.65] 2.60] 1.74] 0.93] 0.83 3.065 2.65) 2.51] 2.31] 1.15} 1.02 
20} 3.99] 2.99] 2.66] 1.83] 1.29] 0.85 5 20 | 3.37] 3.02] 2.64] 2.56] 1.49] 1.20 
3.69} 2.68) 2.51] 1.79) 1.27] 0.78 3.07} 2.69] 2.46) 2.27] 1.45) 1.02 
21); 4.00] 2.90| 2.66} 1.85) 1.33] 0.99 5621] 3.40] 2.92| 2.64] 2.55) 1.58] 1.26 
3.72] 2.64] 2.51] 1.82] 1.30] 0.96 3.13} 2.63] 2.46] 2.27] 1.49] 1.15 
@ 2) 3.96| 2.75] 2.62} 1.96| 1.45] 1.26 5 22 | 3.40] 2.77] 2.66] 2.46) 1.81] 1,41 
3.72| 2.56] 2.48] 1.93] 1.41] 1.26 3.18| 2.54) 2.49] 2.23] 1.67] 1.39 
3.0 | 3.83} 3.00] 2.49] 1.98] 1.65] 1.04 530] 3.26] 3.03] 2.64] 2.30] 1.88] 1.42 
3.56| 2.69} 2.41] 1.98] 1.59] 0.92 3.03] 2.70} 2.48] 2.16] 1.83] 1.16 
3 1 3.84] 2.91) 2.49] 2.03] 1.66] 1.16 5 31] 3.31} 2.91) 2.65] 2.29] 1.94] 1.47 
3.59} 2.63] 2.41] 2.02) 1.61] 1.08 3.09} 2.62] 2.50! 2.16] 1.86} 1.29 
Bp 2} 3.81| 2.71] 2.50) 2:17] 1.68} 1.42 532] 3.31} 2.73) 2.69| 2.26] 2.08] 1.60 
3.60] 2.52] 2.40] 2.14) 1.64] 1.40 3.14] 2.58| 2.47] 2.15| 1.93] 1.55 
3S 3 | 3.69] 2.50] 2.50} 2.40} 1.71] 1.71 533] 3.23) 2.85) 2.50] 2.27] 2.13] 1.78 
3.51] 2.38] 2.38] 2.35] 1.69] 1.69 3.09} 2.73] 2.34] 2.16] 1.98] 1.78 
00} 3.68} 2.99] 2.99] 2.22] 0.82] 0.82 540] 3.05} 3.05| 2.66] 2.26] 2.04] 1.65 
3.23] 2.68} 2.68] 2.09] 0.74] 0.74 2.86} 2.71] 2.56] 2.15| 1.99] 1.29 
10j 3.71} 2.99] 2.87] 2.23} 1.00] 0.88 541) 3.15] 2.86) 2.70] 2.28} 2.03] 1.70 
3.31} 2.68] 2.59] 2.11] 0.96] 0.79 2.95} 2.60] 2.58} 2.16] 1.98] 1.41 
1) 3.73] 2.93] 2.82} 2.24] 1.12] 0.98 5 42] 3.16} 2.83} 2.62] 2.33] 2.01] 1.82 
3.37] 2.65] 2.55] 2.12] 1.05} 0.94 3.01] 2.71] 2.41] 2.19] 1.97] 1.67 
e 0 | 3.70) 3.00] 2.62) 2.25) 1.39] 1.03 543] 3.12) 2.98) 2.38| 2.37} 1.98] 1.97 
| 3.38] 2.69] 2.44] 2.16] 1.37] 0.91 3.00} 2.86] 2.25| 2.21] 1.95} 1.91 
m1) 3.71) 2.91| 2.61] 2.27] 1.46] 1.13 550] 3.07} 2.81] 2.68] 2.50] 1.87] 1.84 
: 3.42| 2.63} 2.43] 2.18] 1.41] 1.07 2.76| 2.73| 2.48} 2.31] 1.84] 1.39 
Cn EE —E—————E———— EEE Ee ee 
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Table 1 (continued) 


4 S 4 ES S S 4 $ 4 S 
PrPyPz| ss Ee Pa ae Pz Py Pz es ry im a 
g 3 3 Hie g 3 ehel tee en 
Digbvak 3.071) 62283!) 122645) 62.50) ol. SE rss 2.80] 2.69| 2.68] 1.76 
286). 2671) O44 | eee oll Soumeeeoe Palos 3.08] 3.03] 2.96} 2.57 
Biome, 3.05| 2.93] 2.50] 2.48] 2.02) 1.81 2.321 2.71 | 2.64), 1.74 
2.95| 2.72] 2.32) 2.30] 1.81 1.76 Tie 0 83.06] 3.04] 2.90] 2.45 
600 3.02} 3.02] 3.00] 2.85] 1.08] 1.08 2S V2.7 Oj S66: eal eT 2 
2.74) 2.68] 2.68| 2.12] 0.92} 0.92 Tao 3.07| 3.02] 2.87) 2.41 
BAL EO SeLOt eS. O02) 9290) ee 7 2S: |e ho DS5)) 27 eek BOs 
9:82) 62-69) | e2i6ai hee. 09) ede OL Om Yes 3.06| 2.97] 2.80] 2.28 
611 S215 lar 2096! | 62:90) e270) ees yo 2.87] 2.69] 2.56] 1.89 
588) | 62.64 li 2. OarieesO 7: \mekecenl Les Pro: 3.06] 3.05 2.73 | 2.22 
Gano 8.09] 3.03] 2.88] 2.53] 1.56] 1.32 2.78) | 2.71) 2158.) 1.98 
2.88| 2.69! 2.60} 2.02] 1.50] 1.09 aoe 3.05| 3.02] 2.70] 2.19 
6.2y1 3.15] 2.92] 2.88) 2.50] 1.68] 1.35 2.83| 2.70| 2.53) 2.03 
9.94) 2.64] 2.59:|%2.01 |) 1.66] 1.21 W302 3.05] 2.94] 2.63] 2.37 
6 Bez 3.14] 2.90] 2.78| 2.38) 1.95] 1.46 2.86| 2.69] 2.43} 2.20 
2.99| 2.62] 2.53] 2.01] 1.73] 1.43 740 3.08] 3.03] 2.48] 2.33 
6 3.0 3.05] 2.98| 2.86] 2.25) 1.94] 1.56 275 te 2s73 52139) 1 2!30 
2.84] 2.71] 2.57] 1.99) 1.84] 1.23 V4] 3.04] 3.02] 2.48] 2.39 
63 1 3.09| 2.91] 2.83} 2.24] 2.03] 1.58 2.81| 2.70] 2.36] 2.34 
2.91) 2.64] 2.56] 2.02) 1.84] 1.35 750 3.09] 3.02] 2.64] 2.16 
6 3.2 3109) W291 e26 72227 eel ole 2.66 9:76 |".2:76 |" *2.52)|) 2:04 
2.97| 2.65] 2.45) 2.14] 1.82] 1.59 8 00 3.13] 3.03] 3.03] 2.59 
Gibko 3.05} 2.93] 2.50] 2.48] 2.02) 1.81 2.81| 2.68] 2.68] 1.62 
DOR MO doi eels OU! MBEcou | ken 810 3.13 | 3.03 |) 2.99] 2.55 
6 4 0 S07 |¢ 2914 9377 (2.31 1.96] 1.81 2.80] 2.69] 2.69] 1.61 
DoT T2718 | Qa eesea ede a) Loe SrleL 3.12 | 3.01 | 2.97 |. 2:50 
641 3.05] 2.91] 2.65] 2.37] 1.95) 1.83 2.80| 2.75] 2.64} 1.59 
2.85! 2.67).2.46| 2.26| 1.77| 1.47 8 20 3.12| 3.04] 2.86] 2.41 
642 3.03] 2.94] 2.52] 2.50) 1.91] 1.90 279) (22370) | P2368 57 
2.92) 2.68] 2.88| 2.29] 1.77) 1.69 821 3.11} 3.01 |) 2.85 |. 2:37 
65 0 3:08) et2 OL i261 ee. 40 2005 7 1 2.78) 2.77) 2.60| 1.64 
2.79) 2.75| 2.42] 2.22) 1.63] 1.46 Sit282 3.10} 2.94] 2.79} 2.24 
66 1 3.05} 2.94] 2.63] 2.387] 2.08) 1.70 Hd le Ow NT Papal elects: 
2.86| 2.69| 2.43] 2.20] 1.63] 1.56 83.0 3.10| 3.06] 2.65} 2.20 
6 6 0 SLO We 2i94-|os2e7 Oe 2i 2a BOT 1.52 27 Or 272 |at2.bS ie Leo) 
2.81) 2.77}. 2.52) 1.88} 1.52) 1:47 SraeL 3.09} 3.02] 2.64] 2.16 
700 3.08]. 3.02:|. 3.02} 2.69) 1.15] 1.15 2776) 2:76:22) 1204 
2.79| 2.68| 2.68] 1.77] 0.96] 0.96 8 40 3:08] 308 |) 12541) 11 52:30 
tend 20) 3.08| 3.03] 3.00} 2.62] 1.29} 1.22 2:74.) 2.74.1) 2540 [192.29 
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Table 2. Values of the frequencies for NaCl. 


For each point the upper value belongs to the approximation b and the lower value to the ap- 
proximation c. py, are the components of the wave vector in units of }- (1/2 ay). 


° 2 Co) 2 3 °° cy ” cy °° co ° 
a 7 T it 7 T T i 7 i ii 1 
—) o So o o i) i) i) —) i] —} o 
‘ 4 4 
ie Pz col col - - - = Py Py De i) = a) o _ ol 
= a *” => 1 2° 4 a ° = 19 2 
3 3 3 3 3 3 3 3 3 3 3 3 


00] 6.02] 3.32] 3.32] 0.00} 0.00] 0.00 422 4.99) 3.17} 3.01) 2.80} 1.97] 1.60 
5.67] 2.95] 2.95] 0.00] 0.00] 0.00 4.69] 2.91] 2.86] 2.63] 1.88] 1.60 

00] 5.95) 3.33] 3.33] 0.71} 0.32] 0.32 430] 4.86] 3.48] 2.85] 2.69] 2.00] 1.77 
5.59| 2.95) 2.95] 0.71] 0.32] 0.32 4.54] 2.99] 2.80} 2.65] 1.88] 1.58 

10]| 5.92] 3.34] 3.27] 0.89] 0.63] 0.46 431 4.87| 3.387] 2.86) 2.77] 2.04] 1.76 
5.58| 2.95] 2.94] 0.89] 0.63] 0.45 4.58| 2.98] 2.81} 2.67] 1.91] 1.66 

y I} 6.89) 3.27] 3.27] 1.06] 0.71] 0.71 432] 4.85] 3.13] 3.04) 2.82] 2.11] 1.82 
5.57| 2.94] 2.94] 1.06] 0.70] 0.70 4.63 | 3.02] 2.81) 2.72] 1.98] 1.79 

00} 5.76} 3.35] 3.35] 1.39] 0.63] 0.63 433 | 4.74] 3.32] 2.91] 2.80] 2.15] 1.96 
5.36] 2.95} 2.95) 1.37] 0.62] 0.62 4.62)" 3.22 | © 2.77" 2.73 |" 2:03) 2.92 

10] 5.75) 3.35) 3.25] 1.45] 0.88] 0.73 440] 4.58] 3.53] 2.88] 2.73] 2.18] 2.07 
5.39| 2.95] 2.93] 1.44] 0.87] 0.71 4.30] 3.02] 2.88] 2.71] 2.03] 1.80 

E 1 5.74] 3.29] 3.23] 1.54] 0.99] 0.88 441 4.60] 3.41] 3.00] 2.73] 2.18] 2.05 
5.41] 2.94) 2.93) 1.52] 0.98] 0.88 4.34] 3.06] 2.84] 2.71] 2.03] 1.88 

2 0 | 5.64) 3.37) 3.11) 1.71)" 1.25) 0.95 are 284.611 93.27 3.15) 25721 - 2:19 | 2.05 
5.382] 2.96) 2.90} 1.71} 1.20] 0.91 4.42] 3.21] 2.78) 2.71] 2.04] 2.00 

21 §.63| 3.29] 3.11] 1.82] 1.28] 1.08 443 | 4.57] 3.53] 2.86] 2.72] 2.19] 2.12 
5.34] 2.93] 2.90] 1.80] 1.24] 1.07 4.49) 3.40) 2.73] 2.71] 2.05) 2.05 

22) 5.54) 3.13] 3.13] 2.10] 1.36] 1.36 454 494.63 | 73.65.) 2.72) 2:72" 2.19) 2.19 
5.29] 2.89] 2.89] 2.08] 1.34] 1.34 4.51] 3.49] 2.71] 2.71] 2.05] 2.05 

00] 5.46] 3.37] 3.37] 2.00} 0.92] 0.92 500] 4.70] 3.44] 3.44] 2.94] 1.39] 1.39 
4.99} 2.96] 2.96} 1.96] 0.88] 0.88 3.94] 2.97} 2.97] 2.80] 1.26] 1.26 

10] 5.46] 3.38} 3.26] 2.04] 1.11] 1.00 510] 4.69] 3.45] 3.35] 2.94] 1.52] 1.48 
5.03| 2.96} 2.93] 1.99] 1.09] 0.96 4.00} 2.97] 2.96] 2.77] 1.44] 1.33 

- I 5.45] 3.31] 3.23} 2.08) 1.24] 1.09 aE I 4.69} 3.38} 3.33} 2.94] 1.71] 1.48 
5.07] 2.94} 2.93} 2.03} 1.21] 1.07 4.05.| 2.98] 2.94] 2.75] 1.60] 1.40 

20] 5.39] 3.40] 3.05] 2.16) 1.48] 1.22 520] 4.64] 3.48] 3.15] 2.91] 1.83] 1.71 
5.05} 2.97} 2.88] 2.13) 1.48] 1.14 4:09)" 2:99) 2:93) 2:71} ° 1.79) 1.52 

21 5.38] 3.31] 3.05] 2.24] 1.54] 1.29 521 | 4.64] 3.38| 3.15] 2.93] 2.01] 1.65 
5.07) 2.94] 2.87) 2.19] 1.48] 1.26 4.13] 2.99} 2.91] 2.70] 1.91] 1.58 

2 2 5.32) 3.14] 3.02) 2.46) 1.67] 1.47 52.2.1 4.61) 3.19 | 3.16) 2.97} 9.24) 1.72 
5.07} 2.88] 2.86] 2.40] 1.61] 1.47 4.21} 3.04] 2.85} 2.68] 2.14] 1.71 

30} 5.18] 3.44] 2.89] 2.40] 1.80] 1.49 653 014.48.) 3.5217 3.07 | 2.77 (92.14) 2.02 
4.89] 2.98] 2.82) 2.40] 1.68] 1.37 4.08) 3.01] 2.96] 2.66] 2.05] 1.75 

3 1 5.18} 3134] 2.89} 2.50) 1.81] 1.54 5 3 1 4.50} 3.41] 3.10] 2.81] 2.25) 1.94 
4.91} 2.95] 2.81} 2.47]: 1.70] 1.47 4.13] 3.06) 2.89] 2.66} 2.12} 1.81 

BS 2} 6.14} 3.11] 2.89| 2.76] 1.85] 1.68 632] 4.51) 3.25) 3.16) 2.87] 2.34] 1.92 
4.93] 2.88] 2.80] 2.68] 1.76] 1.66 4.21] 3.18] 2.80} 2.66] 2.23] 1.89 

33 | 4.99] 3.06] 2.90} 2.90] 1.90] 1.90 533 | 4.46] 3.46] 2.92) 2.90] 2.33] 2.02 
4.85} 2.99] 2.79} 2.79] 1.83] 1.83 4.27) 3.35] 2.74] 2.68] 2.23] 1.99 

Or) 5,09) 3.41)" 3.41) 2.53) 1.17) 1.17 540] 4.23}. 3.68] 3.07] 2.77| 2.34] 2.23 
4.50| 2.96] 2.96] 2.45] 1.09} 1.09 3.92] 3.04] 3.03] 2.65] 2.13] 1.98 

20} 5.09| 3.42] 3.30] 2.54] 1.33] 1.26 541 4.28| 3.45) 3.17| 2.81] 2.31] 2.20 
4.55| 2.97) 2.94] 2.44] 1.28] 1.17 3.98| 3.17] 2.89] 2.66] 2.16} 2.04 

im 1 6.09| 3.35] 3.27] 2.57) 1.49] 1.30 5 4 2 | 4:35] 3.39) 3:17) 2.83) 2.27} 2.16 
: 4.60] 2.95] 2.94] 2.45] 1.43] 1.26 4:10) 3:32) 2:76) 2.67) 2.22) 2.07 
? 0 | 5.03} 3.44] 3.06] 2.60] 1.67] 1.48 543 | 4.42] 3.56] 2.91] 2.81] 2.24] 2.16 
4.61} 2.98] 2.88] 2.49} 1.62) 1.35 4.29) 3.44] 2.71] 2.69| 2.18) 2.06 

i 1 5.03] 3.34] 3.06] 2.65] 1.78] 1.49 550] 3.93] 3.65] 3.08] 2.94] 2.62} 2.11 
4.65) 2.95] 2.88) 2.52] 1.71| 1.44 3.65} 3.07] 3.06] 2.69] 2.18] 2.09 
a RE TS a a ae ee Pe ee 
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Table 2 (continued) 


| oh ab Ae peoaehetes & dovvsp ae ofh eae ap eo giant a 

PyPy Pz) cn cy is = — PzPyPz) “ a ns = a i 
3) 28 b) Mehl edd: ome seen 3 ¢ | 88 01 R824 sega ae 

55 1| 4.05] 3.48] 3.20] 2.93] 2.54| 2.10 3.23| 3.04] 2.99| 2.81] 1.62] Lé 
3.72| 3.23] 2.87| 2.69| 2.25] 2.08 711: 4.11] 3.51| 3.42] 3.23] 1.99] La 

55 2| 4.23| 3.41] 3.17] 2.93] 2.43] 2.07 3.26| 3.09| 2.94] 2.78] 1.81] 1. 
3.93| 3.34| 2.71] 2.70| 2.34] 2.05 720| 4.03| 3.53| 3.40] 3.13] 2.04] 2¢ 

600| 4.37| 3.47] 3.47] 3.20] 1.55] 1.55 3.26| 3.09| 3.01] 2.68] 2.00] 1.1 
3.40| 3.02| 2.98] 2.98| 1.37] 1.37 721) 4.01] 3.53] 3.35] 3.12] 2.31] 1 

610] 4.35| 3.48| 3.41] 3.17] 1.67] 1.65 3.32| 3.11] 2.93| 2.65] 2.17] 1. 
3.47| 3.03| 2.98| 2.90] 1.55] 1.45 722 3.96| 3.50] 3.26] 3.06] 2.62] 1. 

611 | 4.34] 3.43] 3.40] 3.15| 1.88] 1.62 3.38| 3.13| 2.88| 2.59] 2.48] Lé 
3.53| 3.03| 2.94] 2.88] 1.73] 1.51 73 0| 3.87| 3.59] 3.30| 2.93] 2.34] 2. 

620| 4.27| 3.51] 3.32] 3.04] 1.95] 1.89 3.26| 3.12] 3.04] 2.50| 2.32] 14 
3.58| 3.04| 3.00] 2.75| 1.91] 1.64 731 3.89| 3.56] 3.25] 2.95| 2.59] 2.1 

621] 4.27| 3.44] 3.31] 3.04] 2.19] 1.78 3.35| 3.13] 2.96| 2.48| 2.45] 1. 
3.64] 3.06] 2.93| 2.74| 2.06] 1.68 732] 3.93] 3.50] 3.10] 3.09] 2.68] 2. 

622] 4.24] 3.36] 3.21| 3.05] 2.46] 1.81 3.49| 3.18] 2.86] 2.67| 2.39] 2. 
3.73| 3.12| 2.83] 2.67| 2.35] 1.79 740) 3.72| 3.67| 3.18] 2.81] 2.68| 2 

630| 4.11] 3.56] 3.26] 2.84] 2.25] 2.22 3.24| 3.08| 3.08] 2.54] 2.25| 2. 
3.60| 3.07| 3.03] 2.60] 2.12] 1.89 741| 3.82] 3.58] 3.15| 2.97| 2.67] 2. 

631] 4.14| 3.47| 3.26] 2.88] 2.44] 2.05 3.33| 3.16| 2.99| 2.55] 2.28] 2. 
3.66| 3.13| 2.94] 2.59] 2.31] 1.91 75 0| 3.79| 3.62] 3.20] 2.96] 2.67] 2.: 

632] 4.18| 3.35| 3.18] 3.01] 2.51| 2.00 3.24| 3.13| 2.96| 2.50] 2.41] 2. 
3.79| 3.23| 2.80] 2.60] 2.44] 1.97 8 00| 4.14] 3.50] 3.50] 3.31] 1.69] 1.4 

633| 4.23| 3.41| 3.17] 2.93] 2.43| 2.07 3.24| 2.98] 2.98] 2.72] 1.46] 1 
3.93| 3.34] 2.71] 2.70| 2.34] 2.05 810] 4.09] 3.51] 3.47] 3.29] 1.80] 1/ 

640 3.90| 3.63] 3.20] 2.82] 2.55] 2.28 3.23| 3.02| 2.92] 2.69| 1.65] 1. 
3.52| 3.09| 3.06] 2.60| 2.22] 2.12 811) 4.04| 3.53] 3.43| 3.25] 2.02] 1) 

6411 3.99| 3.50| 3.20] 2.96] 2.46] 2.25 3.22| 3.11] 2.94] 2.66| 1.84] 1. 
3.59| 3.20| 2.93| 2.61] 2.27] 2.13 820] 3.96] 3.54] 3.38] 3.19] 2.07] 2. 

642| 4.14] 3.36] 3.22] 3.01] 2.37] 2.19 3.19] 3.10] 3.01] 2.59] 2.03] 1. 
3.79| 3.31| 2.80] 2.64] 2.32] 2.10 8 21| 3.93| 3.56] 3.34] 3.15] 2.36] 1 

650] 3.73| 3.70| 3.12] 2.94] 2.83| 2.06 3.20| 3.18| 2.92] 2.57] 2.21] 1. 
3.40| 3.11] 2.96| 2.70] 2.32] 2.04 8 22] 3.84| 3.60] 3.23| 3.04] 2.70] 1 

651 | 3.93] 3.50] 3.10] 3.09| 2.68] 2.05 3.31| 3.15| 2.82] 2.59] 2.47] 1 
3.49| 3.18| 2.86] 2.67] 2.39| 2.02 8 30| 3.81| 3.60| 3.23] 3.00] 2.37] 2: 

660| 3.84] 3.60] 3.23] 3.04] 2.70] 1.90 3.14| 3.14] 3.05| 2.44] 2.37] 1. 
3.31| 3.15| 2.82] 2.59| 2.47] 1.86 8 31] 3.79] 3.62| 3.20] 2.96] 2.67] 2 

700] 4.19] 3.49] 3.49] 3.30] 1.65] 1.65 3.24| 3.13] 2.96] 2.50] 2.41] 2 
3.22| 2.98] 2.98] 2.88| 1.44] 1.44 8 40| 3.69] 3.69| 3.14] 2.72] 2.72] 2 

7101] 4.14| 350| 3.46] 3.26] 1.77] 1.75 3.09] 3.09| 3.14| 2.51] 2.23] 2 
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